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Abstract. Using Hubble Space Telescope observations with STIS, we study the main exciting stars of N81, a high 
excitation compact Hll region in the Small Magellanic Cloud (SMC). These far UV observations are the first 
spectroscopic measurements of stars in such a region and reveal features characteristic of an 06-08 stellar type. 
The astonishing weakness of their wind profiles and their sub-luminosity (up to ~ 2 mag fainter in Mv than 
the corresponding dwarfs) make these stars a unique stellar population in the Magellanic Clouds. Our analysis 
suggests that they are probably in the Hertzsprung- Russell diagram locus of a particularly young class of massive 
stars, the so-called Vz luminosity class, as they are arriving on the zero age main sequence. 
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1. Introduction 

Understanding the formation of massive stars, which is 
still a largely unsolved problem, requires studying them 
at the earliest phases where they can be reached through 
the enshrouding material at different wavelengths. While 
high resolution radio continuum observations allow the 
investigation of ultracompact H ii re gions formed around 
newborn massive stars (Churchwell 1990| ), high angular 
resolution observations in the ultraviolet, visible, and 
infrared are also necessary to access accurate physical 
parameters of these stars in order to identify their evo- 
lutiona ry stat es (Walborn & Fitzpa trick 1990 , Walborn 
et al. 19951: , Hanson et al. 1996). In particular UV 



observations are of prime importance since massive stars 
emit the bulk of their energy in this wavelength range. 
In practice though observing newborn massive stars is 
not straightforward for several reasons. Mainly, they are 
very rare, and the relatively small evolutionary timescales 
involved make it difficult to catch them just at this 
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very point in their evolution, that is when they become 
observable in the UV and visible (Yorke & Kriigel 1977, 
Shu et al. |l987t Palla & Stabler | 199C| , Beech & Mitalas 



1994 Bernasconi & Maeder |1996| ) 



We have amply argued that the compact H li regions 
known as HEBs (High Excitation Blobs) provide the 
best opportunities for a direct access to massive stars at 
very early stages of their evolution (Heydari-Malayeri et 
al. 2001a| and references therein). The members of this 
distinct and very rare class of ionized nebulae in the 
Magellanic Clouds are small and compact {^-^ 5" to 10" in 
diameter corresponding to ~ 1.5-3.0 pc), in contrast to the 
typical Hii regions in those galaxies, which are extended 
structures (sizes of several arc minutes corresponding to 
more than 50 pc, powered by a large number of exciting 
stars). In general, HEBs are also heavily affected by local 
dust, as one would expect from their very young age 
(Heydari-Malayeri et al. 2001a and references therein, 
see also Israel & Koornneef 1991). And their study is 
pertinent to understanding the process of massive star for- 
mation especially in the context of the Magellanic Clouds. 

Our recent high resolution imaging with the Hubble 
Space Telescope (GO 6563, GO 8246) using the Wide 
Field Planetary Camera (WFPC2) has for the first time 
resolved several HEBs which had appeared featureless 
to ground-based telescopes: SMC N81, N88A, LMC 
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159-5 (the Papillon nebula), N83B, and NllA (Heydari- 
Malayeri et al. 1999a, 1999b, 1999c, 2001a, 2001b). 
The HST observations uncover the so far hidden stellar 
content as well as the nebular features of these compact 
nebulae and display a turbulent environment typical 
of newborn massive star formation sites: outstanding 
emission ridges created by shocks and cavities sculpted 
in the ionized gas by the powerful winds of massive stars, 
prominent dust structures protruding from hot gas. The 
observations also bring to light even more compact H ii 
blobs, immersed in the HEBs, harboring newborn, hot 
stars. 

The present pa per is devoted to N8 1, als o known as 
DEM 138 (Henize |l956t Davies et al. |l976|) , a nebula 
only ~ 10" across and located in the Shapley Wing at 
~1°.2 (~1.2kpc) from the main body of the SMC. A 
first detailed study of this compact Hii region carried 
out by Heydari-Malayeri et al. ( 1988|) , revealed its nature 
and some of its physical characteristics: gas density 
and temperature, chemical composition, mass, age, etc. 
Subsequently, near infrared observations showed the 
presence of H2 emission towards N81 (Israel & Koornneef 



198S), while ^^C0(l-6) emission at two points towards 
this Hii region was also detected (Israel et al. 1993| ). 



However, due to the lack of sufficient spatial resolution, 
it was not possible to view and study the exciting star(s) 
hidden inside the ionized gas. Therefore, the rather 
important question, which is often challenged by star 
formation theories, of whether N81 was powered by 
a single massive star or a cluster of them, remained 
unanswered. This is, however, a critical question for 
theories of star formation. 

High spatial resolution imaging with HST allowed us 
to resolve N81 and reveal the presence of a tight cluster 
of newborn massive stars embe dded in this compact 
nebula (Heydari-Malayeri et al. 1999| , hereafter Paper 
I). Six of the stars are grouped in the core region of 
~ 2" diameter, with a pair of the main exciting stars 
in the very center separated by only 0"27 or 0.08 pc. 
The images also displayed conspicuous marks of strong 
stellar winds, shocks, and ionization fronts characterising 
turbulent massive star forming regions. Moreover they 
revealed prominent dust lanes dividing the nebula into 
three lobes. One of the lanes running over 15" (4.5 pc) 
ends in a magnificent curved plume. A remarkable ab- 
sorption "hole" or dark globule of radius ^ 0"25 {'^ 0.07 
pc) is situated towards the center of the Hii region, 
where the extinction reaches higher values (Ay = 1.3 
mag). These absorption features are probably parts of the 
molecular cloud which has given birth to the massive stars. 



spectral classification for these very young massive stars 
and study their nature. 



2. Observations and reduction 

The General Observer Program No 8246 devoted to 
observations of N81 was performed with Space Telescope 



Imaging Spectrograph, STIS (Woodgate et al. 1998) on 
board HST on 28 and 31 October 1999. The spectra were 
obtained with the far-UV Multi-Anode MicroChannel 
Array (MAMA) detector in the G140L mode covering the 
wavelength range 1120-1715 A. All the observations were 
made through the 52"x 0"2 entrance slit. The effective 
resolution was 0.6A per pixel of 25 //m, corresponding 
to a dispersion of 24Amm~^, or a resolution of 1.2 A 
(FWHM). The exposure times were set according to the 
apparent magnitudes of the stars in order to equalize the 
signal-to-noise ratios (S/N) of the spectrograms. Total 
exposure times varied from 1229 sec (stars #1 and #2) 
to 3169 sec (stars ^^3, #4, and #8). Three relatively 
faint stars (#5, #7, #10), not initially scheduled for 
observations, happened to lie on the slit when observing 
their adjacent stars (#3 and #11). The S/N ratio is 
particularly weak for these stars, yet we present the 
spectrogram of star #5 which shows some interesting 
features. STIS was also used to obtain the spectra of 
the N81 stars in the visible domain. The grating G430L 
covered the range of 2900 to 5700 A with a resolution of 
of 2.73 A per pixel. The CCD pixels of 21 /im yielded 
a dispersion of 130Amm~^. The exposure times ranged 
from 24 sec (star #1) to 750 sec (star #8). 

The calibrated output products from the standard 
pipeline use a default extraction aperture of 22 pixels 
(0"53 on the sky). We carefully reprocessed the 2D images 
using the most recent calibration reference files applicable 
to the observations and extracted the spectra using slits 
of both 6 and 2 pixels. We verified the centering of the 
stars on the slits, and tested for the effects of different 
sky background extraction on our spectra. The 6 pixel slit 
yielded spectra which are very similar to those produced 
by the standard pipeline with an insignificant loss in S/N. 
Even the 2 pixel slit did not indicate any extraction ef- 
fect due to its size other than the expected loss in S/N. 
Comparing the resulting line profiles as a function of the 
aperture width, wc found that other than an increase in 
noise when going to the smaller slits, and slightly modi- 
fying the slope of the spectra (because of slight tilt of the 
spectral images), the line profiles do not change by any 
significant amount. The spectra displayed in Figs. 1 and 
2 are based on 6 pixel extraction apertures. 



From the Stromgren uvby imaging with WFPC2 we 
carried out the photometry of some 50 stars towards 
N81. This allowed us, using color-magnitude diagrams, to 
select the main exciting stars of the region. This paper 
is devoted to the spectroscopy of these stars. We derive 



3. Results 

The N81 stars observed with STIS are listed in Table 1, 
and their physical location can be seen on Fig. 2 of Paper 
I. The table also presents the corresponding photometry 
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Table 1. SMC N81 stars observed with STIS 



Star number 


a 


5 


V 


h~y 


Av 


Mv 




(J 2000) 


(J 2000) 


[y o47M) 














(mag) 


(mag) 


(mag) 


(mag) 


1 


01:09:13.1 


-73:11:38.3 


14.38 


-0.10 


0.22 


-4.84 


2 


01:09:13.0 


-73:11:38.0 


14.87 


-0.11 


0.19 


-4.32 


3 


01:09:13.4 


-73:11:38.4 


16.10 


-0.08 


0.31 


-3.21 


4 


01:09:12.8 


-73:11:38.3 


17.41 


+0.07 


1.02 


-2.61 


5 


01:09:13.3 


-73:11:37.6 


18.29 


-0.05 


0.46 


-1.17 


8 


01:09:12.8 


-73:11:40.2 


17.84 


+0.15 


1.40 


-2.56 


11 


01:09:13.7 


-73:11:33.3 


15.74 


-0.10 


0.22 


-3.48 


13 


01:09:16.1 


-73:11:29.1 


16.65 


-0.08 


0.31 


-2.66 



of the stars (Paper I). The color excesses E{B — V) 
were derived from E{b — y) using the intrinsic color 
(b_—y)o = —0.15 mag for hot stars (Relyea & Kurucz 
1978| ) and assuming that our observed colors represent 
the standard Stromgren system. Then the relation 



E{B -V) ^ lA9E{b - y) (Kaltcheva & Georgiev [1992D 
was used to transform into the Johnson system, which 
finally yielded the extinctions Ay = iAE^B - V). The 
estimated absolute magnitudes are based on a distance 
modulus M-TO = 19.0 mag (corresponding to a distance 
of 63.2 kpc, e.g. Di Benedetto 1997 and references therein) 
and assuming that the Stromgren y filter is equal to the 
Johnson V . 



The final reduced spectrograms are presented in Figs. 
1 and 2, where the former figure includes the four bright- 
est stars of the sample, whereas Fig. 2 displays the fainter 
ones. The main stellar features (Cm A 1176, NvA1239, 
1243, OvA1371, SiivAA1394, 1403, Civ A 1548, 1551, 
and He II A 1640) are distinguished with tick marks. The 
labels appearing below the features indicate cases where 
contamination with an interstellar component is possible. 
The identification of S v A 1502 A is based on the work of 
Werner & Ranch ( ^OOlD . 



An outstanding aspect of these spectra is the extreme 
weakness of the UV wind profiles. Weak stellar wind fea- 
tures in SMC O stars have already been found by sev- 
eral workers (Hatchings 1982, Garmany & Conti 1985, 
Walborn et al. |l995b ), who ascribed them to the metal 
deficiency of the SMC leading to a reduced radiation pres- 
sure responsible for driving the winds of early type stars. 
More recent observations have further confirmed this re- 



sult (Smith Neubig & Bruhweiler 1997, Walborn et al. 
2000| ). However, the wind features observed in the stars of 
N81 are even weaker. If we consider the usually stronger 
wind lines seen in O stars (such as the N v A 1239, 1243 and 
C IV A 1548, 1551 in dwarfs, or the Si IV AA 1394, 1403 in gi- 
ants/supergiants), with the exception perhaps of star #5, 
none shows the emission part of a P-Cygni profile and the 
absorption is extremely weak particularly for N v AA 1239, 
1243. For stars which we classify as O types (see below), 
such a behavior is observed for the first time. 



3.1. Sp ectra I classifica tion 

Traditionally the spectral classification schemes of stars 
are based on their optical spectra. However, recent studies 
comparing spectral features in the optical and UV have 
resulted in a global, co heren t picture of classification 
criteria (W alborn et a l. 1985 , 1995a, Smith Neubig 
Bruhweiler 1997t 



199£). In particular. Smith Neubig 



Bruhweiler (1997) have proposed a UV classification 
system for O and B stars of the SMC which is defined 
by a set of standard, low resolution spectra observed 
with the International Ultraviolet Explorer (lUE). This 
UV scheme, which was used by the authors to derive 
classifications for 133 O and B stars of the SMC, while 
independent of the MK system, shows general agreement 
with those deduced from visual data. 

The low S/N ratio of our STIS optical spectra and 
in particular their contamination with strong nebular 
emission lines limit their practical use for spectral 
classification. This can be understood since N81 is a 
very bright compact H ii region with strong nebular 
emission lines in the visible part of the spectrum (Paper 
I). Although nebular emission lines are present also in 
the UV part, they are much less troublesome. Therefore, 
we will use the method put forward by Smith Neubig & 
Bruhweiler (1997). However, given some morphological 
differences of the N81 UV spectra with previously well 
studied stars, the limitations of the optical part of 
the spectrum mentioned earlier, and the constraints 
available from the UV classification scheme (Smith 
Neubig & Bruhweiler 1997), it is clear that there is no 
unique solution to the spectral classification of our targets. 

The presence of He ii A 1640 and O v A 1371 in all 
spectra, except perhaps in stars #13 and #8 which have 
lower S/N ratios, is the first and strongest evidence that 
we are observing O type stars. These spectra though 
display characteristics suggesting that the stars also 
belong to the dwarf luminosity class since the features 
Si IV AA 1394, 1403 as well as NvAA1239, 1243 and 
CivAA1548, 1551 are weak (even weaker than usual). 
Star #5 may be a different case as explained below. These 
features are known to increase with luminosity, ranging 




Fig. 1. Rectified HST/STIS ultraviolet spectrograms of the four brightest stars in SMC N81. The prominent absorption 
feature at A 1210 A is due to the Lya.The wind profiles arc indicated with tick marks and the features possibly 
contaminated by an interstellar component are labelled below the lines. 
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Fig. 2. Spectrograms for the remaining four les bright stars of N81. The notation is the same as in Fig. 1. 
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from weak P-Cyg profiles on the main sequence to very 
pronounced P-Cyg profiles in the supergiants (Walborn 
et al. [19954 Smith Neubig & Bruhweiler |1997D . A dwarf 
luminosity class is also supported by the optical spectra 
which show no He ii A 4686 and N ill A 4640 emission. 
Although the morphology of the N81 spectra differs 
qualitatively from that of the known O types, we may 
classify them as "zero age main sequence dwarf" Vz, 
based on the weakness of the wind lines (Walborn & 
Parker 119921 Walborn & Blades |1997l Walborn et al. 



200C). Note that the original definition of the Vz class is 
He II A 4686 absorption much stronger than He ii A 4541 
or He I A 4471, and therefore the association of these stars 
with that class is indirect. 

No clear distinction between the various subtypes can 
be made based on the available spectral features, but it is 
very Hkely that all stars are of a late O type (~ 06-08). 
This is supported by the following three facts. First, the 
presence of the OvA1371 feature indicates a spectral 
type earlier than ^ 08 (Smith Neubig & Bruhweiler 
1997). Second, the weakness of the Sim A 1300 feature, 
which appears in the wing of Oi + Siii A 1304, excludes 
much later types. Finally, the N v AA 1239, 1243 feature is 
weaker than CivAA1548, 1551, which is only seen in O 
types > 06-07 or alternatively in OC stars (Walborn & 
Panek p^ ). 



Star #5 shows some puzzling emission features in its 
spectrum (AA 1480, 1508, 1616, and 1640 A), at least one 
of which. He ii A 1640, is apparently part of a P-Cyg pro- 
file. It should be stressed though that the S/N ratio is not 
high enough to be absolutely certain about their presence. 
If these were indeed wind induced features, their presence 
would confirm the fact that the CivAA1548, 1551 pro- 
file has a particularly marked emission component com- 
pared to the other stars of the sample. Since both of the 
He II A 1604 and S v A 1502 fines show P- Cyg p rofiles in Of 
and Wolf-Rayet spectra (Walborn et al. 1985, Wilfis et al. 
1986| ), star #5 appears to be an Of or WR candidate in 



their structure in detail. 

Using an effective temperature derived from the 
estimated spectral t ypes (06-08), in particular the one 
from Vacca et al. ( |1996 ) for dwarfs, and the absolute 
magnitudes from Table 1, we place the N81 stars (filled 
circles) in an HR-like diagram (Fig. ||). For further 
comparison we also include the mean My magnitudes 
for 03-BO dwarfs based on a compilation by Vacca et 
al. (1996), as well as the mean My relation of Walborn 



( |1972D used by Walborn & Blades ( |1997[ ) for 30 Doradus 
stars. Fig. ^ clearly shows that most of the stars in N81 
are sub-luminous compared to the mean My for dwarfs. 
This finding still holds even if the spectral types are 
shifted by 1-2 subtypes towards later types. Compared to 
the so-called Vz luminosity class, the sub-luminosity of 
the N81 stars is more pronounced. Our targets are up to 
~ 2 mag fainter in My than the mean relation for dwarfs! 
Although this is not the defining characteristic of the Vz 
class, the above indications as well the weakness of the 
UV wind lines further attests that these stars belong to 
the Vz class. 



Two sets of ZAMS and isochroncs of 4 Myr are pre- 
sented in Fig. H using the Geneva stellar evolution tracks 
for metallicities bracketing approximately that of the SMC 



N81. 



(Lejeune & Schaerer 2001). One set has been computed for 
a metallicity 1/20 Solar {Z = 0.001) and is marked with 
a long-dashed line, while the second has a metallicity of 
1/5 Solar {Z = 0.004) and is marked with a dotted line. 
The observed My and spectral types are roughly com- 
patible with positions close to the theoretical ZAMS or 
young ages. Given mostly the lack of an accurate spectral 
subtype determination we cannot firmly establish if the 
lower luminosity stars (#8, #13) are really hotter than 
the ZAMS. Atmospheric modeling is in progress to obtain 
more accurate stellar parameters of these unique young 
stars in the SMC (Martins et al., in preparation). Based 
on the Z — 0.004 tracks, the ZAMS luminosities and 
masses corresponding to the observed My are between 
logL/L© ~ 4.2-5.5 and ~ 14-50 Mq respectively. 



3.2. Stellar parameters and wind properties 



4. Discussion 



To constrain the stellar wind properties we have examined 
the line profiles of the strongest UV lines. The best indi- 
cations for velocity shifts come from the C iv A 1548, 1551 
feature. The profiles of stars #1 and #2 show blue-shifted 
absorption reaching up to velocities of ~ 1700 and 2000 
km s~^ respectively. These values are compatible with 
terminal velocity measurements in other SMC O stars 
of similar spectral type (Walborn et al. 1995 
& Crowther 1998). For star #5 we derive a 
velocity of 1000 km s~^, while in the other profiles the 
terminal velocities are lower. However, these profiles also 
show asymmetries which are stronger on the red side, and 
further investigations will be necessary to understand 



Prinja 
terminal 



The HST spectra of N81 presented here are the first ones 
ever obtained from a tight cluster of stars in a HEB. 
The reason is that these stars, embedded in a compact 
emission nebula, have not been reachable by ground- 
based telescopes. And even with recent developments in 
ground-based instruments, taking spectra of individual 
stars in the visible remains still practically infeasible. As 
a result, contrary to other massive stars in the SMC that 
ha ve been obser ved from space in the UV (Walborn et 
al. |l995bi ^OOOj Smith Neubig & Bruhweiler |l997| ), the 
N81 stars lack high quality spectra in the visible. Our 
low-resolution HST spectra in the visible, imposed by 
stringent time allocation constraints, were intended to be 
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4.7 4.6 4.5 



Fig. 3. The absolute My magnitude versus effective temperature diagram of the N81 stars is compared to the zero 
age main sequence (ZAMS) and 4Myr isochrones at vario us metaUicities. The TcS of the stars has been estimated 
using an 06 or 08 spectral type and the Vacca et al. (1996) scale. The ZAMS and 4 Myr isochrones for a metallicity 
of 1/20 Zq (Z = 0.001) are indicated by long-dashed lines while the same pair for 1/5 Z© (Z = 0.004) is plotted with 
dotted lines. The ZAMS curves are not plotted for log (Tcff) > 4.7 due to lack of an appropriate conversion to My- 
Also shown, as solid lines, are the My- Toff calibration for dwarfs from Vacca et al. (1996) as well as the one from 
Walborn (1972). Note that four of the observed stars are situated in the HRD-like diagram in a locus suggesting that 
they are either on the ZAMS or that they have a young age. The lower luminosity stars (# 4, 8, 13) are possibly 
hotter than the ZAMS. 



a first exploratory such attempt. 



Massive stars observed from the ground and also 
with HST are typically much brighter than the ones seen 
in N81 and comparably they are much less affected by 
nebular emission and dust. The decoupling from nebulae 
is presumably due to the evolutionary state of these stars; 
they have had enough time to entirely and/or effectively 
disrupt their Hii regions and the associated dust. This 
means that those bright stars are older than the N81 



members, as supported by their spectra. Among the 15 



SMC stars studied by Walborn et al. (2000), 9 are clearly 



giants, 5 are peculiar and have already developed emission 
line features of NivA4058, N ill A 4640, or Cm A 4650 
(two of them being on the main sequence), and the last 
one is a pure main sequence. 

The fact that all the observed exciting stars of N81 
display the Vz characteristics, further supports the very 
young age of the cluster (Paper I). This observed concen- 
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tration of Vz stars in a small region is intriguing since in 
the LMC region of 30 Doradus only 6 of the 104 O and 
early B stars classified by Walborn & Blades (1997), that 
is ~6%, belonged to the Vz category. If we assume that 
all 20 bright stars we detected towards N81 (Paper I) are 
of O or B type, then the lower limit for the fraction of 
Vz stars in N81 is nearly 35%, which is quite considerable! 

The exact nature and evolutionary stage of Vz stars 
is still unknown. Presumably these objects arc close to a 
transition from their formation locus in the HR diagram 
to the main sequence. However, several issues regarding 
their properties, which are relevant to our observations of 
compact low-metallicity H ii regions, remain open: 



1. Why do these stars show such weak stellar winds? Are 
their mass loss rates compatible with expectations 
from "normal" O stars, i.e. due to their reduced 
luminosity compared to stars with similar effective 
temperatures and representative for massive stars in 
their earliest evolutionary stage? 

2. Are these objects truly on the ZAMS, blueward or 
redward of it? The presence or absence of massive 
stars on the ZAMS - corresponding to the locus of 
completely homogeneous objects initiating H-burning 
- yields information on the star formation process. 
Indeed, the apparent lack of Galactic OB stars close 
the ZAMS could be due to the hiding of such stars 

198^ ) 



in their parental cocoon (Garmany et al. 
explained by the progressive redward bending of the 
upper part of the birthline due to moderate mass 
accretion rates in an accretion scenario for these stars 



(Bernasconi & Maeder |1996| ). In the former case the 
position of the "earliest" star visible provides thus 
constraints on the duration of the hidden phase. For 
the latter scenario, the position of the bluest stars 
constrains the accretion rate Mace; the existence of 
massive stars close to the Z AMS requires high values 
of M ace (Norberg & Maeder |2000i Behrend & Maeder 
|2001| ). 

Given the present uncertainties on the Tcff determi- 
nation of the N81 stars, one may speculate that some 
of the objects are indeed hotter than the ZAMS, as 
indicated by Fig. ^. If true, this "blue straggler" Hke 
behavior cou ld be indicative of stellar collisions (e.g. 
Benz & Hills 1987), and thus a signature of formation 
of massive stars via this process, as advocated among 
others by Bonnell et al. (|998|). 



3. Could some of the N81 stars still be accreting pre- 
main sequence objects? The stars with smaller My 
magnitudes (# 3, 4, 8, 11, and 13) have estimated 
luminosities slightly above the probable high lumi- 
nosity Herbig Ae/Be pre-main sequence stars found 
by Lamers et al. (1999) and de Wit et al. (2001) 



in the LMC and SMC. Their location in the HRD 
also coincides with the region where the predicted 
birthlines follow quite closely the ZAMS. Given the 
strong indications for a very young age of N81, it is 
thus conceivable that our objects are still accreting 
mass as part of their formation process. The redshifted 
C IV profiles in the above cited objects could be an 
indication of ongoing accretion. 

4. Are Vz stars related to low metallicity? Vz stars have 
also been detected in the Milky Way by Walborn and 
co-workers (private communication). Therefore, it ap- 
pears unlikely that they are due to a simple metallicity 
effect. 

While at present these answers remain fairly specu- 
lative, our upcoming determinations of stellar and wind 
parameters and further studies of similar objects should 
hopefully shed more light onto these issues related to the 
formation and early evolution of massive stars. 

5. Conclusions 

Our HST/ STIS observations of the brightest massive 
stars powering the high excitation compact Hii region 
SMC N81 reveal that the stars have strikingly weak wind 
profiles and a pronounced sub-luminosity which are clear 
indications of their early evolutionary state. Most likely 
they belong to the Vz category of massive Magellanic 
Cloud stars (Walborn & Parker 1992| , Walborn & Blades 
1997), which are very young stars just arriving on the 
ZAMS or already located near it. These stars may 
also serve as templates for newborn massive stars of 
distant metal-poor galaxies which cannot be individually 
resolved. Therefore, a more detailed study and modeling 
of their properties is highly desirable as it will help shed 
some more light on the intricacies and consequences of 
the early stages in stellar evolution. 
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